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Abstract 

For the linear sigma model with quarks we derive renormalization 
group flow equations for finite temperature and finite baryon density 
using the heat kernel cutoff. At zero temperature we evolve the ef- 
fective potential to the Fermi momentum and compare the solutions 
of the full evolution equation with those in the mean field approxi- 
mation. We find a first order phase transition either from a massive 
constituent quark phase to a mixed phase, where both massive and 
massless quarks are present, or from a metastable constituent quark 
phase at low density to a stable massless quark phase at high density. 
In the latter solution, the formation of droplets of massless quarks is 
realized even at low density. 
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1 Introduction 



Recently a transparent approach to the evolution of a hadronic system with 

resolution has been constructed using a Schwingcr proper time representation 
of the fluctuation determinants of quarks and chiral mesons [1]. This method 
gives a picture of the transformation of constituent quarks at low resolution 
into partonic massless quarks at high resolution observed in deep inelastic 
scattering [2] . At finite temperature the same type of renormalization group 
(RG) flow equations give critical indices for the chiral phase transition in 
agreement with the 0(4) model [1, 3]. Present finite temperature QCD 
lattice simulations seem to indicate such an 0(4) type behavior, with some 
uncertainty. 

Simulation of QCD at finite baryon density on a lattice is still a challenge: 
In Euclidean space the chemical potential gives rise to a complex action which 
forbids Monte Carlo calculations. Quenched simulations at finite density suf- 
fer from additional shortcomings [4]. The challenge is therefore to come up 
with a calculational scheme for finite baryon density which has a well con- 
trolled predictive power. Theoretical studies of high density matter given in 
this work are indispensable to understand what is going on with high-energy 
heavy ion collisions which probe not only high excitation energy (perhaps 
temperature) but also high baryonic density. Recent studies have suggested 
a very rich phase structure at high baryon density [5, 6]. There has been 
considerable work devoted to extract the equation of state of nuclear matter 
in terms of nucleon - nucleon potentials. The problem is to link the high den- 
sity region accessed by high energy heavy ion collisions with the low density 
nuclear physics region. 

At low density it is probably not very efficient to describe baryonic matter 
by quarks, but in an intermediate region, quarks and mesonic bound states 
may be the right degrees of freedom. The mesons are formed by the strong 
gluonic attraction in certain channels. The chiral linear sigma model with 
quarks is a good model to investigate the dynamics below a momentum scale 
of A; < kuv — 1.2 GeV^. Such a hybrid approach has been already used at 
finite temperature with success. In this paper we will treat pure quark matter 
and concentrate on chiral symmetry aspects of the transition. We neglect 
quark confinement but assume the gluons to be confined into the mesons. 

^The definition of this momentum scale depends on the renormahzation group scheme. 
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The connection to nuclear matter will be considered elsewhere [7]. The finite 
baryon density dynamics is much more sensitive to neglecting confinement 
than the finite temperature dynamics. Indeed it is known that e.g. the 
NJL model gives a first order phase transition at unrealistically low baryon 
density [8, 9]. Without confining quark forces, which repel quarks tripled in 
neutral nucleons from each other [10], the linear sigma model with quarks 
probably overestimates the binding energy per baryon of nuclear matter. 

We use renormalization group flow equations in the presence of flnite 
quark density. As a first step, we start with the heat kernel representa- 
tion of the effective potential using the cut-off function of the earlier finite 
temperature calculation [1]. The extension of this technique to finite den- 
sity is straightforward in contrary to the cut-off function introduced by the 
Wetterich group [11] for finite temperature studies. However, the proposed 
method has its shortcomings: 1. The heat kernel cutoff is a Lorentz invari- 
ant cutoff, and as such does not zoom to the Fermi surface with increasing 
evolution steps. In principle the goal of the renormalization group approach 
for finite fermion density is [12] to treat the critical long wave length particle 
hole excitations of the finite Fermi system towards the end of the evolution. 
These particle hole excitations are generated at the Fermi surface. In the 
relativistic case with Goldstone bosons this final evolution step should co- 
incide with the treatment of the zero mass bosons. 2. Extra terms which 
reduce the ultraviolet sensitivity in our heat kernel cut-off function result in 
an integration pole at the Fermi surface, making the numerical integration 
difficult. 

In this paper we concentrate on deriving the formalism for a finite den- 
sity calculation with finite temperature. We solve the fiow equations for finite 
density at temperature T = leaving our final goal, the determination of the 
phase diagram in the whole (T, /x) plane to another paper. Furthermore, we 
compare our renormalization group results with the mean field approxima- 
tion. 

As a description of choice in this work we use the linear sigma model. We 
think that this model can play the role of a simplified theory showing the 
effects of relativistic field theory at finite density as a kind of Ising model 
for nuclear physics. The equation of state improved by the renormalization 
group in the one loop gives clear signals on what is happening in the finite 
density system. The main message is the intricate connection between the 
effective mass of the fermions and the effective masses of the mesons ex- 
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changed between the fermions. In classical nonrelativistic nuclear physics 
the time scale of the meson exchange is short in comparison with the time 
scale of the nucleon dynamics. At normal nuclear density the mesons will not 
be modified drastically by the nuclear medium. In the quark picture the time 
scale for the exchange of a meson is not different from the time scale con- 
nected to the rearrangement of the fermion momenta. The mesons cannot be 
taken into account as potentials, which are turned into effective interactions, 
due to the quark rescattering in the medium. The mesons participate fully 
in the reorganization of the chiral symmetric phase with increasing density. 
The quark many body problem is therefore very different from the nuclear 
physics problem where the nucleon dynamics largely decouples form the me- 
son dynamics. The only possible tool we have to solve such an intrinsically 
nonlinear problem in field theory is the renormalization group, therefore it is 
worth to study this technique in nuclear physics to get a better understanding 
of the phase structure of baryonic matter. 

The paper is organized as follows: in section 2 we review the derivation 
of the renormalization group flow equation for the effective potential at zero 
density and temperature. The extension of this method to finite temperature 
and baryon density is presented in section 3. At the end of this section we 
show a first result from the RG flow equation at finite density. In section 4 we 
discuss the mean field approximation to the linear sigma model and in section 
5 we present our results for two sets of mean field couplings, obtained from 
the RG evolution. Here we study the finite baryon density phase transition 
of the linear sigma model in the mean field approximation. In section 6 
we compare the results of the grid and the mean field calculations. Section 
7 is devoted to a summary and to the conclusions. In the Appendices we 
discuss the connection between the fermionic part of the fiow-equations and 
the mean field result and give a detailed derivation of the fiow equations at 
finite temperature and density. 

2 Evolution of the linear a model 

We consider the chiral constituent quark model with quarks, a and tt mesons. 
At zero temperature and chemical potential in Euclidean space the partition 
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function or the generating functional without external sources A, is given by 
Z[A = 0] = J VqVqDaV7rex.p{- J d'^x (Cp + JCb)} (1) 
with a fermionic Cp, and a bosonic jCb, parts, 

Cj, = q{x) YieOe + g{(7 + ir7f75)] q{x) , (2) 



Cb = \[{d,af + {d,T:f]+F{a^ + i:^). (3) 

The Yukawa coupling of the constituent quarks to the mesons is denoted by 
g. The parameters of the linear cr-model at T = are chosen in the same 
way as in ref. [1, 3]. We assume that at an ultraviolet scale kuv, the full 
QCD dynamics reduces to a hybrid description in terms of quarks and chiral 
bound states. The gluons are assumed to be frozen in the residual mesonic 
degrees of freedom and their couplings. 

At the beginning of the evolution at k^v = 1-2 GeV we choose the effective 
potential density Fq, to be of the following form: 

fS) = '^<t>' + \<l>\ with </>2 = a2+7^^ (4) 

where the positive mass squared m? = 0.4^ GeV^ reflects a symmetric ground 
state, i.e. the minimum of the potential lies at the origin. The four boson 
couphng at this scale is A = 30 [1, 3]. The values correspond to a critical 
temperature ~ 150 MeV, and a chiral symmetry breaking scale k^r ~ 1 
GeV which one obtains after performing the RG procedure. 

The effective potential density F{(j)), can be evolved [1] using the heat 
kernel method. For this purpose the one loop effective potential is calculated 
with a cutoff function /(/c^r), which contains the evolution scale k: 

f{x = ^t) = e-^ (l+x+^x^^ . (5) 

Doing this, the couplings of the of the effective potential become scale de- 
pendent: m = m{k) and A = \{k). 
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Using the the Schwingcr proper time representation the fermionic 
contribution to the effective potential density is: 



^ ' 





and correspondingly [1] the bosonic part is: 

2J T ■'^ ' J (27r^ 
^ ' 



1Y e-'^[9^+ff^</''^l 



(6) 



(27r)^ 



(7) 



with 



{2ny 



Fo^Fo{<l>\k), F^ 



902 



d'F, 



(8) 



The total effective potential density is the sum of these two terms, F — 
F^ + F^. The evolution equation of this potential results from the derivative 
of the potential with respect to k. In the spirit of the renormalization group 
improved 1-loop approximation the derivative only acts on the cutoff function 
and the potential density Fq, is replaced by the evolving potential density 
F. The evolution equation for the hnear sigma model then has the following 
simple form: 



dF 
'dk 



k^ 



1 



327r2 [k^ + 2F' fcH2F'+4F"02 ^+^2^2 



(9) 



In the approximation to the RG-evolution used here only the effective po- 
tential density F, evolves with the scale k. The Yukawa- coupling g = 3.23 
is assumed to be constant. The limitation of this approach will be discussed 
later. Note that during the evolution one passes from the region with > 0, 
where the potential is symmetric, to the region < where the potential 
has a mexican hat shape. The denominators of the meson loop terms (e.g. 
the first two terms on the l.h.s. of eq. (9)) indicate a limitation of the one 
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loop renormalization flow. The one loop corrections are of order 0{h) in 
an expansion in h. In the regime k'^ + 2F' < however, the usual Gaussian 
fixed point is unstable and a non-Gaussian fixed point should be considered 
bringing in 0{h^) effects [13]. With our choice of parameters this happens 
eitk< kinst ~ 200 MeV. 



3 Evolution equation at finite density 

The renormalization group flow equations give a well determined shape of 
the effective potential with which we will work in this section. Since they 
are formulated in the continuum, the finite baryon density pb can be imple- 
mented. Although the primary aim of the present work is to discuss the chiral 
transition at finite baryonic density, we shall present the master formulae for 
finite temperature T and finite pb and then take the limit T — > 0. 

We start our derivation from eqs. (6,7) extending them to finite temper- 
ature and finite chemical potential. Since we are working now with fixed 
temperature and chemical potential, one should replace the potential density 
F, by the thermodynamical potential density Q{(f)'^). The fermionic part of 
the effective potential density generalized to finite temperature T and chem- 
ical potential p is 

= ^AN^Nj fik^r) E /-^e-[(-+^'^)^+^-^+^^^% (10) 

with the cutoff function f{x) given by eq. (5), and the Matsubara frequencies 
f„ = (2n + l)7rT for fermions. The bosonic part is the one already examined 
in [1], 
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with Matsubara frequencies cun — 2mrT for bosons. The derivatives in the 
potential are taken again with respect to 0^. The total effective potential 
density is the sum of the two terms, Q = + . 

An advantage of the heat kernel regulator is that the calculation of deriva- 
tives with respect to the momentum scale k can be performed analytically 
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to yield compact formulae. As shown in appendix A, we have 



Stt- 



d 



^k'-- I dq 



E, 



q,k 



(12) 



where Eq^k = \/q'^ + k'^ + g'^(f)^ and n{x) and n{x) are the Fermi-Dirac dis- 
tribution functions for particles and anti-particles. Evaluating the derivative 
with respect to k"^ on the right hand side, we obtain 



327r2 
1 



r dq 



1 



Elk 



[1 - n{Eg^k) - n{Eg^k)) 



{n(E,,fe)(l - n(E,,,)) + n(E,,fe)(l - n(E,,,))} 



(13) 



We notice that the distribution functions explicitly show how the temperature 
and the baryonic density modify the T = and = result: The right 
hand side of the evolution equations in vacuum is diminished (1) by the Pauli 
blocking effect as seen in the first bracket, and (2) by the thermally excited 
states as seen in the second bracket. Similarly, the bosonic part at T 7^ is 
evahiatcd as shown in appendix A, where dO,^ /dk^ is expressed in terms of 
the Bose-Einstein distribution function. 

Although formula (13) is generally valid for T 7^ and 7^ 0, it is 
not useful when the temperature is set to zero due to the factor T^^ in 
the second bracket. The zero temperature formula is obtained easily from 
eq. (12) making the substitutions n{Eq^k) — > 9{ijl — Eq^k) and n{Eq^k) — > 0: 



dk 



NcNf 



k^ 



87r2 A;2+^202 



\/ ix^ — k'^ — g^(f)- 



^e(//-7^^T^) 



(14) 



Here again we see how the evolution equation at finite baryon density is 
modified by Pauli blocking. In appendix B we show how the integral of the 
finite density part of equation (14) gives the mean field result, which will be 
presented in section 4. 

For the charge neutral system (zero bosonic chemical potential) at T = 
the bosonic term is the same as in vacuum. Adding both contributions, we 
have the zero temperature flow equation for finite baryon density as follows: 



'dk 



k^ 



1 



327r2 [k'^ + 2n' A;2 + 2fi'+4nV 



(15) 
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We stress that such a compact form is obtained for the flow equation in 
the heat-kernel method; the Pauli blocking effect due to the presence of the 
Fermi sea is explicitly represented by a theta-function. This is certainly an 
advantage of the method. However, the presence of the theta-function makes 
numerical evaluations difficult because its derivatives produce singular terms 
which are not easy to control in the numerical analysis. 

As an attempt to circumvent this difficulty, we let the chemical potential 
run during the evolution; so the chemical potential will be a function of the 
scalar field and the momentum scale /i — ii{(f), k), as proposed by Shankar in 
reference [14]. In fact, it turns out that no singularities of the fermion-terms 
appear in the evolution with a running fi. In the functional integral for Q we 
explicitly insert a /c-dependent chemical potential. 

The method is best explained starting from k — 0. At the end of the evo- 
lution the transition to the free energy density can be made via a a Legendre 
transformation. 

F(A; = O,02) = Q(A; = O,02)+p//(A; = O,02), (16) 
where /x(fc = 0) has to be eliminated from fl{k = 0, 0^) via the equation 

dn{k = 0) 



d^{k = 0) 



-p. (17) 



Now we make a small change of infrared cut-off scale k and adjust /x(/c) in 
such a way that the density remains constant: 

-p. (18) 



dp{k) dp{k) 

The same procedure is repeated at each step of k. The explicit evaluation of 
the l.h.s. at a finite evolution scale k using (10) yields 

2Ar / 3,2 3 

with 



xf = Jp2 - A;2 - ^202_ (20) 
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An analysis of the solutions jJi^xp) from eq. (19) shows that at /c = this 



equation has a unique solution (i = \J g'^^P' + k'p^ with the Fermi momentum 
kp defined by p = 2Nckp/37r'^. For A; 7^ it is advantageous to discuss the 
solutions in terms oi z — xp/k. Equation (19) then reads 

P-^^k'-9iz), (21) 

with 

^(^) = ^' + ^^ + ^- (22) 
The behavior of g{z) is shown in Figure 1; g{z) has the minimum 

_ 1 a/3 + 1 

'\/3- 1 



at 2; = = Zm- At finite evolution scale k the line of constant 

•?^2„ .„ . , . ...... 



2N, 



will in general cut the constraint function at two points Zi and Z2- It 



is necessary for the equation to have real solutions that 

1/3 



..l^^l ,24) 
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At A; = kth the line of constant j^-js is tangential to g{z) ai z = z„i. 
Corresponding to the two solutions, the chemical potential has two different 
values at fixed density: 



j^^l^k) = ^Jk^Zi{ky + k^ + g^(l)^ (^ = 1,2). (25) 

The derivatives d^i/dk are calculated by demanding that the baryon density 
found at /c = remains the same independently of the evolution scale k. Thus 
the k derivative of the l.h.s of eq. (19) is zero, leading to 

d/ii 1 k 
'dk ~ 'J^iSzf + 4:zf -1' 

Applying the chain rule one can calculate the flow equation for the thermo- 
dynamic potential with running chemical potential. 



(26) 



'dk 



on 

'dk 



+ 



ti=const. 



dn dfi{k) 
dfj, dk 



dQ 
'dk 



-p- 



^=const. 



dii{k) 
' 'd^ 



(27) 



or rearranging this formula for each k we get the equivalent free energy at 
each k: 



dF 
~dk 



d{n + pfi) 



Pb 



dk 



dn 

'dk 



(28) 



fj,=const. 



We can eliminate the chemical potential in the flow-equation at each k via a 
5-f unction. 



/ df 5 [f{z)] , 



with 



2N, 



''k'g{z)-p. 



37r2 



(29) 



(30) 



The evolution equation (27) with running chemical potential has two terms 
i — 1,2 corresponding to the two roots of the constraint equation: 



dn 

'dk 



dfS[f{z)] 
df 



=1,2 •' 



dz 



dz 



dn 

'dk 
df 



dz 



+ 

fi=const. 
-1 

S{z-Zi] 



dn dn{k) 
djjL dk 



dn 

'dk 



IJ,=const. 



dpjk) 

'W 



(31) 
.(32) 
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Figure 2: Normalized energy per baryon plotted against inverse density 
, normalized to normal nuclear density po = 0.16 fm""^, in the bro- 
ken phase (dashed line) and in the symmetric phase (solid line). The 
evolution eq. (33) is solved on the grid. The energy of the broken phase 
at zero density is eo. 



Note that at the two zeros of the function /, Zi and Z2, the Jacobians have 
opposite signs (see figure 1) leading to some cancellation. Thus we arrive at 
the final flow equation 

^ ^ I ^ ] (33) 



dk 



The reason of the presence of the theta-function in (33) is apparent from (10). 
Note that the term l/(8zf+4:zf—l) is singular at k — ktu but can be integrated 
analytically. 

Equation (33) is the evolution equation with running chemical potential 
nik) at zero temperature for the thermodynamic potential density. In order 
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grid calculation 




1 2 3 4 5 

Pb/po 



Figure 3: The binding-energy per baryon calculated from the evolution 
equation on the grid is shown as a function of normalized to normal 
nuclear density po- 

to follow the minimum of the order parameter at fixed density we need the 
free energy density. Recall that at A; = there is only one solution to the 
constraint equation (19) which is related to Z'z, therefore we can calculate the 
free energy density F aX, k — Q from the thermodynamic potential density 
Q,{k — 0) and the chemical potential ^{k — 0): 

F(A; = 0,(/)^) = Q(A; = 0,(/)^) +p//(A; = O,0^) (34) 

Due to the theta term in eq. (33) the k evolution of Vl{k) and F{k) is 
the same till kth and does not feel the baryon density. Below kth the density 
effects set in with the contributions of a shell of fermions. During the course 
of evolution the outside radius of the shell increases and the inside radius 
diminishes until the longest wavelengths in the fermi sea arc integrated. The 
effect arising from the fermi sea is classical, hence 0{h^) and dominates the 
RG flow equation. For large enough densities p 0.45po, we have kth > hnst 
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and eq. (9) is replaced by the density driven evolution without unstable boson 
terms. 

We solved eq. (33) numerically on a grid: the full potential density is 
discretized as a function of 0^ on a grid of hundred points between < 0^ < 
0.05 GeV^. The resulting hundred differential equations are solved with a 
Runge Kutta method. 

As we discussed earlier at the end of section 2 the meson terms in the 
flow equation (33) develop singularities due to the instabihty of the effective 
potential. This behavior is well known in the literature [15, 16], the mesonic 
effective potential in the one loop approximation generates tachyonic meson 
masses. The renormalization group scale k"^ avoids these poles for some part 
of the evolution, but cannot get rid of them totally. The singularities in the 
boson denominators of the evolution equation appear at kf^g^-\-2Q! — 0, indi- 
cating the disappearance of the Gaussian fix point [13] . As our studies show 
the inclusion of the non-Gaussian fixed point does not change the position 
of the minimum considerably at zero temperature. Hence in this paper we 
neglect the singular mesonic contributions beyond k = kinst in the evolu- 
tion equation leaving the more exact solution to a future work. The results 
of the grid calculation are shown in figures 2 and 3. We will discuss them 
extensively in section 6 and compare them with the mean field calculations 
presented in the next section. 

4 Coarse grained potential in mean field ap- 
proximation 

Since the numerical solution of the RG fiow equations has difficulties at small 
momentum scale k, we discuss another approximation to the low momentum 
region in this section. We evolve the vacuum theory from the ultraviolet scale 
k — kuv down to A; = A;^,' corresponding to normal nuclear matter density 
Pb = Po- Thereby wc obtain a coarse grained potential which is appropriate 
for the dynamics at these low momenta. This potential contains the vacuum 
loop effects of the quarks and bosons integrated out up to this scale. Then 
we solve this coarse grained linear sigma model in mean field approximation. 

To keep the mean field approximation transparent we approximate the 
coarse grained potential with the original fourth order form of equation (4) 
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Figure 4: The potential 1^(0^) at momentum scale k = kth, from the grid 
calculation (solid line) and the fitted curves with mean field parameters 
X (dashes) and Y (dashed-dotted line). See table 1. 



with renormalized parameters. These parameters depend on the range in 0^ 
where the fit to the coarse grained potential is done. We have chosen two 
sets: a fit on a wider region in 0^ (X) and a narrower one (Y). The resulting 
parameters are summarized in Table 1. The coarse grained potential and 
the two fits are shown in figure 4. The negative values of the mass squared 
indicate that at the scale k = kp '^e already entered the broken phase, with a 
nonvanishing expectation value of the a field, a. A mean field solution with 
these potentials is straightforward by evaluating the Hamiltonian in the 
presence of the Fermi sea and then minimizing H with respect to a. Since we 
have already taken into account quantum fluctuations from the Dirac sea and 
the bosons in the evolution, we do not need to consider these pure quantum 
fluctuations any longer. Their effects are assumed to be of higher scale than 
kp and therefore they are integrated up in the coarse grained potential. Of 
course, we are missing typical quantum many body fluctuations not included 
in the mean field approximation of the many body system. 
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Type 


Fit range 0^ 


2 








SV4 




X 


[0,0.05] 


-0.260 


30.0 


0.0940 


0.902 


0.15554 


0.720 


Y 


[0,0.01] 


-0.082 


9.18 


0.0945 


0.906 


0.11633 


0.405 


grid 








0.0944 


0.906 







Table 1: Effective meson potentials used. The fit range and the mass pa- 
rameter m\^p are in GeV^, the ctq, 5^^^ and the mass of the a meson are in 
GeV. 



The mean field solution of the resulting 0^ theory is standard. The Hamil- 
tonian reads 



H 



I 



(fx 



{dp + gpa) q + 



2 \ 

— — H : — u 



(35) 



Using a plane-wave basis for the quarks one can rewrite the Hamiltonian as 
an integral over momenta. In the ground-state, T = 0, the Fermi-sphere is 
filled from the bottom up to the Fermi-momentum kp- The energy density 
is evaluated to be 



E 



4A^c 



kp 



MF -4 



(36) 



47r2 



4-4 1 

-g a log 



^/j4 + g^^ + kF\ 1 



A, 



go- 



+ 2^mf(^ + 4 



-a . 



The quark Fermi-momentum is fixed by the quark-density p = 3pB, 

kp - 



2N, 



c 



(37) 



The mean cr-field configuration is calculated by minimizing (36) with 
respect to a. One ends up with a self-consistent equation for the mean field 
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There is always the trivial solution cr = corresponding to the symmetric 



phase. At zero density the non trivial solution = y— m^/A has a lower 
energy and represents the spontaneous broken phase. At higher density, the 
chirally symmetric phase has lower energy. The energy density in the broken 
phase relative to the symmetric one is 

4 

60 = SbripB = 0) = -— = -B, (39) 

whereby we define S as a kind of bag constant. It gives the amount of 
energy density by which the partonic vacuum lies above the constituent quark 
vacuum. We estimate the sensitivity of the calculation to the input coarse 
grained potential by comparing the minima of the energy per baryon in both 
phases. The massless partonic (symmetric) phase has an energy density 

^.^-imr^r (40) 



with a minimum of the normalized energy per baryon 



37r \ / m \ 



2iVe / I A 



(41) 



PB 

at 

In the constituent quark phase the asymptotic value of the energy per baryon 
is three times the quark mass: 

^hr{p) — ^0 



PB 



= 3^ao . (43) 



Therefore in order to have a stable broken phase at low density the following 
condition is necessary (but not sufficient), 
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Wc note, that with g=3.23 this condition is fulfilled for parameterization {X) 
but not for (Y). As a result the latter does not have a stable homogeneous 
phase with broken chiral symmetry. 

In the following section we present our numerical results. 



5 Results in mean field approximation 

In this section we discuss the results of the mean field calculations, shown 
in the figures 5-7, before we compare them with the grid calculation in the 
next section. In figure 5 we present the energy per baryon as a function of 
the inverse density relative to the energy of the broken phase at zero baryon 
density. In the low density limit, the energy per baryon of the broken phase 
slowly approaches three times the constituent quark mass Sttlq = SgaQ. This 
is the preferred state at low densities, where we have a noninteracting dilute 
system of constituent quarks. Since our coupling constant g, is fixed from the 
beginning to yield a constituent quark mass of 300 MeV, the vacuum mass of 
the "nucleon" is smaller than 938 MeV in both cases (X) and (Y) (cf. table 
1). In the chiral limit which we pursue here the nucleon is lighter than the 
real one. In the upper left corner of both plots the parabola represents the 
partonic phase. 

For the couphngs (X) condition (44) is fulfilled, thus the phase diagram 
may be obtained using the Maxwell construction (dashed line): 

For a first order phase transition the equilibrium-condition for a given 
temperature (in our case T = 0) and pressure is 

G(T, P, N) = min . (45) 

In the region of coexistence between the two phases (I/partonic) and (Il/con- 
stituent quark) , the temperature, pressure and chemical potential have to be 
equal to each other. Hence 

fi'/' = -^G{T, P, Nb) -^iF + Py)-^ + (46) 

Nb JylB J^B PB 

where the pressure P is constant in the coexistence region. At zero tem- 
perature the free energy per baryon F/Nb equals the energy per baryon 
E /N — e/ pb and we can read off from the tangent construction the energy 
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mean field calculation (X) mean field calculation (Y) 




Figure 5: The energy per baryon plotted against in units of nor- 
mal nuclear density = 0.16 fmT^ for mean field parameterization X 

(left) and Y (right). Note the different energy scales in the figures. The 
dashed line in the left plot represents the Maxwell-construction which 
determines the region of coexistence between broken and chiral symmet- 
ric phase. 



per baryon as a function of the inverse baryon density in the mixed phase: 

Nb Pb 

The slope is the negative pressure and the intercept with the vertical axis 
is the baryon chemical potential at the phase transition. The dashed line 
in figure 5 connects the low density constituent phase with the high density 
partonic one. The phase transition takes place between 0.27 and 1.90 times 
normal nuclear density. In between the two phases coexist. 

The couplings (Y) do not fulfill condition (44), e.g. we have no stable 
broken phase. Hence there is no Maxwell construction in this case. 

For both parameterizations (X) and (Y) we calculate the pressure from 
equation (47) . The result is shown in figure 6. It shows clearly the stable 
broken phase for parameterization (X). For parameterization (Y) there is 
no stable broken phase. At low densities the pressure is always negative, 
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Figure 6: The pressure as a function of inverse baryon density, calcu- 
lated for mean field parameterizations X (left) and Y (right). The thick 
dashed line in the left plot is the Maxwell-construction. In the case Y 
(right) no Maxwell-construction exists. Note that the negative values of 
the pressure are solutions of the equations but have now physical inter- 
pretation. In both cases the pressure for the symmetric phase (a = 0) 
asymptotically approaches the bag pressure B (cf. table 1) for small 
densities (dashed-dotted line). 



thus it is more advantageous to pack the quarks into droplets of massless 
quarks with nonperturbative vacuum between them than to have loosely 
bound constituent quarks associated with nuclcons. 

In figure 7 we present the binding energy per baryon subtracting the 
"baryon mass" 3mQ, from the energy per baryon. Note that in the case (X) 
the system is unbound, the energy of the partonic phase has a minimum at 
~ 1.8 times normal nuclear density. This minimum lies on the edge of the 
coexistence region. In the case of couplings (Y) the system is strongly bound, 
the minimum of the partonic phase is Pb = 0.78 po- 

The difference between the two cases {X) and {Y) demonstrates that 
the form of the coarse grained potential Q infiuences the physics at low mo- 
mentum scale drastically. In the mean field approximation it is mainly the 
sigma mass = 2\aQ which determines the amount of attraction i.e. bind- 
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mean field calculation (X) mean field calculation (Y) 




1 23 4 50 1 2 3 4 5 

Pb/ Po Pb/ Po 

Figure 7: Binding energy per baryon is plotted against in units of 
normal nuclear density po, for the couplings X (left panel) and Y (right 
panel). 



ing or unbinding of baryonic matter. Recall the attraction in nonrelativistic 

2 2 

Hartree approximation varies as — f;^- Potential (X) with a high cr-mass of 
0.728 GeV is less attractive than the potential (Y) with a u-mass of 0.404 
GeV. The meson-meson interaction term A, also fixes the structure of the 
intermediate density region. A large A gives a mixed phase as produced by 
potential (X) cf. eq. (44). The grid calculation shares the low field strength 
region with the potential ( Y) , and has the mass of the a meson in between of 
the ones from the potentials (X) and ( Y) (see figure 4) . Therefore we expect 
that the equation of state on the grid lies between the extremes determined 
by the potentials (X) and (Y). 



6 Comparison of mean field results with grid 
calculations 

In this section we compare the results obtained from the mean field approxi- 
mations with the ones from the grid calculation. The main results are shown 
in figures 8 and 9. 
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Figure 8: Normalized energy per baryon, calculated on the grid (dashed 
line) and in the mean field calculations with couplings X (upper solid 
curve) and Y (lower solid curve), is shown as a function of inverse density 
Pg^, normalized to normal nuclear density. 



In figure 8 we present the energy per baryon obtained in the (X) param- 
eterization of the mean field (upper curve), the (Y) parameterization (lower 
curve) and the grid calculations (dashed line) . The three curves show a very 
similar behavior. At low density the energy per particle is approaching the 
same limit of three times the mass of the constituent quarks, as we discussed 
in the previous section. The nuclcon mass is ~ 900 MeV in both cases. With 
increasing density, one arrives at the point where no broken phase is sup- 
ported any longer, i.e. the mean field eq. (38) does not have a nontrivial a 
solution and only the massless partonic phase exists. For the coarse grained 
couplings (Y), this happens at ps = 0.32 poi while for the grid calculation at 
Pb = 0.56 Pq. The parameterization (X) leads to highest transition point at 
Pb = 1.32 po- In all three cases the order parameter drops to zero, i.e. we 
have a first order phase transition. 

Since the minimum of the energy per particle in the partonic phase is 
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Figure 9: Binding energy per baryon, calculated on the grid and in 
the mean field approximation with potentials X (upper curve) and Y 
(lower curve), is shown as a function of normalized to normal nuclear 
density. 



below the minimum of the broken phase, no Maxwell-construction is possible 
and the broken phase is not stable even at lower densities for both the grid 
calculation and model (Y). This phase starts from zero density and persists 
until Pb = 1.17po foi' the grid and until ps = 0.78 po for the mean field 
calculation (Y). At these densities the droplets fill up the whole volume. They 
are bound by ^ 74MeV (grid) and ^ 172 MeV (mean field), respectively, per 
baryon number. If ones compresses the system above this density the fermi 
pressure pushes the equation of state higher up in energy. 

Contrary to these cases there exists a Maxwell construction for potential 
(X) as we have seen in the previous section. The region of coexistence be- 
tween the stable constituent quark phase and the partonic phase ranges from 
0.27po to 1.90 po- 

In figure 9 we show the equation of state for the three calculations. The 
coarse-grained couplings (Y) as well as the grid-calculation yield bound sys- 
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tcms. However, in both cases the binding is too strong, 74 MeV per baryon 
in the grid and 172 MeV per baryon in the mean field calculation with model 
(Y). The saturation densities lie near normal nuclear density, but baryonic 
matter is already in a chirally symmetric phase. 

7 Summary and Conclusion 

We have calculated a coarse grained effective potential from renormalization 
group fiow equations in a quark model with explicit meson fields. At the 
fermi momentum scale we continue the evolution including Pauli blocking. To 
avoid meson instabilities we switched off the meson loop terms at kfh ^ kp. 
The resulting quark matter overbinds and is in a chiral symmetric phase. 
We presented two mean field approximations to the full grid calculation with 
different sets of couplings. The grid solution of the flow equation, with the 
evolution of the meson effective potential omitted below kth, lies in between 
the result of these two mean field approximations. 

The many body physics found here is very similar to the NJL model where 
one finds a first order phase transition with either a mixed phase or a droplet 
phase [9]. Also the instanton induced quark interaction of ref. [5] produces a 
droplet phase of partonic quarks. In the effective potential the determining 
equation is the relation between the energy per quark in the partonic phase 
and the constituent quark mass (cf. eqs. (41,43)). 

The shape of the droplets in the low density phase have to be determined 
from an independent calculation including surface effects. The bag produced 
by the linear sigma model does not confine. It has a finite height, outside the 
massless quarks acquire a constituent quark mass. Otherwise the solution of 
the evolution equation is similar to a MIT bag type solution with massless 
quarks inside. Up to now attempts to model the nucleon as a soliton in the 
linear sigma model have failed due to the instability of the sigma solution 
arising from the integration over the sea quarks. The evolution equation may 
be helpful in finding such a soliton solution too. 

In our approximation we used lowest order in two cases to address the fi- 
nite density problem: the coupling constant g, was fixed during the evolution 
and the wave function renormalization Z was set to one. A recent publica- 
tion of the Wetterich group [11] uses a running meson-quark coupling and 
obtains a solution with a mixed phase extending from very low to very high 
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density. In fact the evolution of the couphng constant at finite density is not 
necessarily the same as at zero density. The equations for the wave function 
renormalization and the coupling are modified at finite density. Important 
effects such as pion condensation [17] may appear after the wave function 
renormalization. One therefore has to make more extended calculations to 
know the full result for the linear sigma model at finite density. 

In quark matter the linear sigma model with only attractive sigma mesons 
and pions is not sufficient to capture the nuclear physics. Repulsive effects 
and confinement into nucleons play an important role. To monitor the transi- 
tion to the deconfined phase one needs an order parameter which keeps track 
of this transition. Since the linear sigma model with quarks is already a 
hybrid model, it is not unnatural to include also nucleon degrees of freedom 
explicitly and trace the transition of nucleon to quark degrees of freedom 
directly. This has been done in a separate paper [7] . 
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A Calculation of dVt^/dk and dVt^/dk 

In this appendix, we calculate the derivative of the thermodynamical poten- 
tials in the heat kernel method. 

The fermionic part of the flow equations: 

The fermion part of the thermodynamical potential is written as 

Vl^ = NcNfI^{k), (48) 
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with 



(49) 



Using the fact that 



n,q 



dfjk'r) 
dk 



(50) 



and putting E'^^^, = g^cjp' + (f + k"^, the derivative dl^ /dk is evaluated to be 

7-1 



dk 



■n,q 



with 



+ 



(51) 



(52) 



where use of the fact has been made in the last equality that the sum w.r.t. 
the Matsubara frequencies does not change when the sign of the frequen- 
cies is changed. 

Now the nice point is that the sum in S{q, k) can be performed analyti- 
cally, as follows; 



S''{q,k) 



d 
dx^ 

d 
dx^ 



E 



ln[(i/„ + ijjif + x^] + ln[(i/„ - ijjif + 



^ ln[i/2 + (a; - nf] + ln[i/2 + [x + fif 



Then evaluating the derivative, we have 

Eq^k — 



S'{q,k) 



E. 



-1 oo 

E 



'<l,k n= 



+ 



pl + {Eq,k-^)^ iyl + {E,,k + fi) 



(53) 



(54) 
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Now utilizing the formula 

oo 

E 

r 

we end up with 



n=—oo n 



1 _ 1 X 



(55) 



q,k — A* 



tanh — '^'^„ ^ + tanh 



Eq^k + 



2T 



'^TEg^k I 2T 
Inserting S^{q,k) into (51), we have 

where n{x) {n{x)) is the Fermi-Dirac distribution function, respectively 

1 . . 1 



(56) 



1 - n{Eg^k) - n{Eg^k) 



(57) 



n(x) — 



, n{x) 



q{x+ij)/T _^ ]_ • 



(58) 



Now let us calculate the derivative w.r.t. /c^, which may be put into the 
integral; since k"^ appears only in the combination q'^ + k"^, the derivative can 
be converted to the one w.r.t. q^. Then making a partial integration, we 
obtain 



k' 



87r2 dk^ 



dq 



1 



E. 



,k 



:i-n{Eg,k)-n{Eg,k)). 



(59) 



Inserting (59) into (48), we finally reach the formula presented in the text; 
T = ^Z^l^^^ I dq- 



dk'^ 



(60) 



The bosonic part of the flow equations: 



Similarly, with the fermion part, the boson part of the thermodynamical 
potential involves the integral; 



n'{k)^-^j:J^ifik'r)e 



-Tlujl+f+m?] 



(61) 



n,q 



27 



with Un being the Matsubara frequencies for bosons. Here we shall con- 
fine ourselves to the case where the boson has no chemical potential. The 
extention to the case with finite chemical potential is easy. 

The derivative w.r.t. k can be performed as much the same way as the 
fermion part. The only difference comes in with the formula 



E / , =-=rcoth4;- (62) 



Thus we obtain, 



dq-—coth—^, (63) 



with = q'^ + k"^ + m'^ . 

Performing the derivative as before, we have 

dfl^ k^ 



dk^ 647r2 



(64) 



+ 2 dqi^nB{Eg,k) + Y^nB{Eq^k){l + nB{Eq,i 
where nB{x) is the Bose- Einstein distribution function given by 



)) 



Taking the limit T — > 0+, we obtain 

dn^ k^ 



dk'^ 647r2m2 + A;2' 



(66) 



B Connection between the mean field and 
the flow-eq. result 

Let us regard only the density dependent part of the evolution equations 
(eq. 14). If this part decouples from the meson evolution, we can integrate 
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this equation analytically: 



''initial W"' 

^ ' (68) 



StT^ i/*2-fl202 + ^202j ^^.^2 _ /j2 _ ^202 ' 

The integral gives the following result: 

+2(sV)Mog^i+^^^Z?!). (69) 
Since one can do the replacements, 



the potential takes the following form: 

AT, /„. —3 



= ^(2fc^^V4+(?202 _ ^2^2 ^^^^2 +^2^2 



4,4, y5±ZZ±^^ AT, 8 / 3 

Recalhng the relation between Q and F 

O = F - 3pB/x , (71) 

we can identify the term in brackets with F, and this is indeed the fermionic 

part of our mean field result eq. (36). 

We can also do this independent check for the equation (33), but because 
of the complicated form of our constraint equation this can not be done 
analytically. Numerical integration of the density dependent part in (33) 
reproduces again the mean field results. 
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